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Abstract. We present HST images at 622 nm and 300 nm of the jet in 3C 273 and determine the run of the optical 
spectral index at 0'.'2 along the jet. The smoothness of spectral index changes shows that the physical conditions 
are varying smoothly across the jet. There is no correlation between the optical flux and spectral index, as would 
be expected for relativistic electrons suffering strong cooling due to synchrotron emission. We find no evidence 
for localized acceleration or loss sites. This suggests that the spectral shape is not changing much throughout the 
jet. We show that relativistic beaming and/or sub-equipartition magnetic fields cannot remove the discrepancy 
between light-travel time along the jet and the lifetime of electrons emitting optical synchrotron radiation. We 
consider this further evidence in favour of a distributed electron acceleration process. 
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1. Introduction 

While radio jets are a common feature of radio galaxies 
and quasars, optical emission has to date been observed 
from only about 15 extragalactic jets. As shown by polari- 
metric observations (starting with Baade (1956 ) for M87), 
both t: 
tinuun 



le ladiu and uplical einissiuii is syiichiuLiun cun- 



quired for the emission of high-energy (optical and UV) 
synchrotron radiation have a very short lifetime which is 
much less than the light-travel time down the jet body in, 
e. g., 3C 273. Observations of optical synchrotron emission 
from such jets (Roser & Meisenheimer 1991, 1999) as well 
as from the "filament" near Pictor A's hot spot (Roser & 



lauiaUun. VVliile inluiiiiaLiuii cm Llie suuiee s inag- 

netic field structure may be obtained from the polarisation 
structure, the diagnostic tool for the radiating particles is 
a study of the synchrotron continuum over as broad a 
range of frequencies as possible, i. e., from radio to UV or 
even X-ray wavelengths, and with sufficient resolution to 
discern morphological details. 

The radio and optical emission observed from hot spots 
in radio jets can be well explained by first-order Fermi ac- 
celeration at a strong shock in the jet (the bow shock) 



(Meisenheimer & Heavens 1986 
19871; iMciscnheimcr et al. 1989| , 



Heavens & Meisenheimer 



1997). But it is not clear 



that the optical synchrotron emission from the jet body, 
extending over Lens of kiloparsecs in some cases, can be 
equally well explained by acceleration at strong shocks in- 
side the jet. As is well known from standard synchrotron 
theory, electrons with the highly relativistic energies re- 
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* Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract 
No. NAS5-26555. 



Meisenheimer 1987; Perley et al. 1997) suggest that both 
an extended, "jet-like" and a localized, "shock-like" accel- 
eration process are at work in these objects in general and 
3C273's jet in particular (Meisenheimer et al. 1997). The 



extended mechanism may also be at work in the lobes of 
radio galaxies, where the observed maximum particle en- 
ergies are above the values implied by the losses within 
the hot spots ( Meisenheimer 1996) and by th e dynamical 
ages of the lobes (Blundcll fc Rawlings 2000). 



The fundamental question is thus: how can we explain 
high-frequency synchrotron emission far from obvious ac- 
celeration sites in extragalactic jets? Although most of the 
known optical jets are very small and faint ( Scarpa & Urry 
20001 ), there are a few jets with sufficient angular size and 



surface brightness to be studied in detail: those in M87 (a 
radio galaxy), PKS 0521 — 365 (an elliptical galaxy with 
a BLLac core), and 3C273 (a quasar). 

We have embarked on a detailed study of the jet in 
3C 273 using broad-band observations at various wave- 
lengths obtained with today's best observatories in terms 
of resolution: the VLA (in combination with MERLIN 
data at A6 cm) and the HST. Using these observations, we 
will derive spatially resolved (at 0'.'2) synchrotron spectra 
for the jet. 3C273's radio jet extends continuously from 
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Filter 


Mean A 


FWHM 


Exposure 


Point 


source 


Extended 


source 




nm 


nm 


s 


mag 


^Jy 


mag/D" 


fiJy/n" 


F300W 


301 


77 


35 500 


26.1 


0.04 


20.8 


5.2 


F622W 


620 


92 


10 000 


27.7 


0.04 


23.0 


2.8 



Table 1. Passbands and limiting magnitudes for the observations. Point source: 10cr detection limit, prediction by 
the exposure time calculator. Extended source: 5tr per pixel detection limit, determined from the background noise 
measured on reduced frames. Magnitudes are Vegamags referred to the corresponding HST filter band. 



the quasar out to a terminal hot spot at 21'.' 5 from the 
core, while optical emission has been observed only from 
10" outwards|] On ground-based images, the optical jet 
appears to consist of a series of bright knots with fainter 
emission connecting them. So far, synchrotron spectra 
have been derived for the hot spot and the brightest knots 
using ground-based imaging i n the radio (Conway et al 
1993|)7|near- infrared if '-band (Neumann et al. 1997 ) and 



optical I, R, B-bands ( Roser fc Mciscnhcimcr 1991) at a 
common resoluti on of 1'.'3 ( Meisenheimer et al. 1996a ; 
Roser et al. 2000). The radio-to-optical continuum can be 
explained by a single power-law electron population lead 



( Biretta et al. 2000| ) on board the Hubble Space 
Telescope^ in three sets of observations on March 23rd 
and June 5th/6th, 1995, for 35 500 s through filters F300W 
(ultra-violet, U) and 10 000 s F622W (roughly R c in the 
Kron-Cousins system). The resulting limiting magnitudes 
are listed in Table |[ The jet was imaged onto the center 
of the Planetary Camera (PC) chip which has 800 by 800 
pixels with a nominal size of 0'.'0455 projected on the sky. 
The telescope was oriented such that the position angle 
of the PC chip's y-axis was « 154°. This way, the quasar 
itself is also mapped on the PC, while each of three neigh- 
bouring astrometric reference stars (see Table 3 in Rosei 



a high-energy cutoff frequency decreasing from 10 17 Hz to 
10 15 Hz outwards along the jet. The aim of the study is 
both the determination of the spectral shape of the syn- 
chrotron emission, and by fitting synchrotron spectra ac- 



ing to i constant radio spectral index 2 of —0.8, but with & Meisenheimer 1991) was observed on one of the three 



cording to Meisenheimer et al. (1989), deriving the maxi- 
mum particle energy everywhere in the jet. The observed 
spectra can then be compared to predictions from theo- 
retical work. 

As an intermediate result of our study, we present HST 
WFPC2 images of the jet in 3C 273 in a red (F622W) 
and near-UV (F300W) broadband filter. The near-UV ob- 
servations constitute the highest-frequency detection of 
synchrotron emission from 3C 273 so far. (Extended X- 
ray emission has also been observed, but it is unclear at 
present whether this, too, is due to synchrotron radiation 



flRoser et al. 2000| ; [Marshall et al. 200l| ; |Sambruna et al 
200l|).)| From these images, we construct an optical spec- 
tral index map at 0'.'2 resolution. 

After a description of the observations and the data re- 
duction in Sections || and ||, we examine the direct images 
in Sect. ||. The creation and description of the spectral in- 
dex map follow in Sect. |[ We analyse the map in Sect. | 
and conclude in Sect. [?]. Details regarding the alignment 
of HST images are found in App. |A|. 



2. Observations 

Observations were made using the Planetary Camera of 
the second Wide Field and Planetary Camera (WFPC2) 



For the conversion of angular to physical scales, we as- 
sume a flat cosmology with f2 m = 0.3 and Ho = hao x 
60kms _1 Mpc~ 1 , leading to a scale of S^/ijjj/kpc per second 
of arc at 3C273's redshift of 0.158. 
2 We define the spectral index a such that /„ oc v a . 



Wide Field chips. This choice gives distortion-free images 
of the jet, while allowing the use of the quasar and the 
reference stars for image alignment. 

The total exposure time was split into individual expo- 
sures of about 2500 s or one HST orbit, executed in three 
separate sets of observations. Each of these exposures was 
done at a telescope pointing slightly offset from all others 
(by up to 1") to facilitate the correction for chip artifacts. 
One short exposure was obtained in each set of observa- 
tions and filter to measure the positions of the quasar and 
reference stars, most of which are saturated on the long 
exposures. Our study of 3C 273's jet will be conducted cm- 
ploying a 0'.'2 effective beam size. This results in a typical 
signal to noise ratio (S/N) per resolution element in the 
red band of around 100 in the brightest regions and 30-40 
in inter-knot regions. Because both the jet flux and the 
WFPC2 and telescope throughput decrease towards the 
UV, the UV-band S/N is only about 40 in the brightest 
regions, while the inter-knot regions are barely detected. 

3. Data reduction 

The data were re-calibrated under IRAF with the STSDAS 
package provided by STScI, using the "most recent" cali- 
bration reference files as of December 1998, because some 
of the files had changed since the original observations. 
The values of background noise measured on the cali- 
brated frames agree well with the values expected from 
photon statistics, as calculated from the expected read 
noise, dark current and sky background level. 



1995 



One of the PC chip's charge traps (Whitmorc & Wiggs 
Voit et al. 1998|) lies inside the jet image, in column 



339. This has no observable effect on the faint UV image, 
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but the effect had to be corrected on the well-exposed 
red-band images. This was done by replacing the affected 
portion of each image by the corresponding pixels from 
an offset image, as a correction according to Whitmore & 



Wiggs (1995) unduly increased the noise in the corrected 
part of the image. 

The images were initially registered using the com- 
manded offsets to the nearest pixels. This alignment is 
sufficient for the rejection of cosmic rays as these only af- 



The criss-cross pattern visible in regions CI and C2, 
and less clearly in Bl-2 and Dl, is reminiscent of a (dou- 
ble?) helical structure ( Bahcall ct al. 1995|), but could also 
be explained by oblique double shocks ( jHardce fc Norman 



1989) 



The jet has three "extensions" (Fig.|l|), none of which 
has been detected at radio wavelengths. The morphol- 
ogy of the outer extension supports the classification 
as a galaxy based on its colours made by Roser & 



feet a iimall number of adjacent pixels. Cosmic rays were 



rejected using a standard K-cr algorithm, rejecting all pixel 
values deviating more than 4cr from the local (low-biased) 
median in a first pass, and neighbouring pixels with more 
than 2.5cr deviation in a second pass. The number of pix- 
els treated this way agrees with the expected cosmic ray 
hit rate for the images. 

A model of the sky background and "horizontal smear" 
(increased pixel values in rows containing saturated pixels 
from the quasar's core, Chap. 4 of |Biretta et al. 2000| ) was 
fitted in the part of the image containing the jet using 
second-order polynomials along rows. The coefficients of 
the polynomials were then smoothed in the perpendicular 
direction. The conversion from photon count rate to phys- 
ical flux units used the throughput information provided 
by STScI which is valid after the 1997 SYNPHOT update. 

4. Maps of the optical brightness 

The calibrated images are presented in Figs, [l] and |[ The 
morphology of the jet is identical in both images and ap- 
pears rather similar to the morphology in high-resolution 



radio maps (Conway et al. 1993; Bahcall et al. 1995). The 
exception to this is the radio hot spot, being the dom- 
inant part in the radio but fairly faint at high frequen- 
cies. Our images show structural details of the optical jet 
which were not discernible on earlier, shallower and un- 



dersamplcd HST WF images of O'.'l pixel size ( Bahcall 
et al. 1995). Based on our new maps, the term "knots 



seems inappropriate for the brightness enhancements in- 
side the jet, as these regions are resolved into filaments. 
The higher resolution necessitates a new nomenclature for 



the jet features (Fig. fl 


) . For consistency with earlier work 


(Lelievre et al. 198' 


1: 


Flatters & Conway 1985; Roser & 


Meisenheimcr 1991' 


, our nomenclature is partly at vari- 



(1995) 



ance with that introduced by Bahcall et al. 

The jet is extremely well collimated - region A has an 
extent (width at half the maximum intensity) of no more 
than 0'.'8 perpendicular to the average jet position angle of 
~ 222° (opening angle ^5°). Even region D2/H3 is only 
1" wide (opening angle w 2.5°). The optical jet appears 
to narrow towards the hot spot, in the transition from H3 
to H2. 

Region A is now seen to extend further towards the 
core than previously known. It may be noteworthy that 



Lelievre et al. (1984 ) reported the detection of an exten- 
sion of knot A towards the quasar, whose existence at the 
reported flux level was not, however, confirmed by later 
work. 



Meisenheimcr (1991). The nature of the other two exten- 
sions, however, remains unknown even with these deeper, 
higher resolution images. The northern inner extension 
was already resolved into two knots (Inl, In2) o n a Faint 
Object Camera image ( Thomson et al. 1993 ). The two 
knots are extended sources and clearly connected to each 
other. The southern extension is featureless and an ex- 
tended source. 

Comparing the direct images, we can immediately es- 
timate that the jet's colour slowly turns redder outwards 
from region A. The similarity of the jet images in both fil- 
ters shows that there are no abrupt colour changes within 
the jet. For a quantitative assessment of the extensions' 
and the jet's colour in the following, we derive an optical 
spectral index map. 

5. Optical spectral index map 

Firstly, we consider how to derive a map of the spectral in- 
dex ctRu from the presented images, at a common effective 
beam size of 0'.'2. 

5.1. Definition of optical spectral index 

Independently of whether a spectrum actually does fol- 
low a power law over any range of frequencies, a local 
two-point spectral index can be defined between any two 
surface brightness measurements B\ , B2 at frequencies V\ 
and V2 (with V2 < v\), respectively, as 



ln^-' 



(1) 



The error on the spectral index is computed from the 
values of the noise ci, 02 in the respective input images: 



1 



In ^ V B f 



(2) 



This formula shows that the spectral index error depends 
on the S/N of the input images and the "baseline" be- 
tween the wavelengths at which the observations are made. 
In the following, we will use v\ = v(U) for the UV and 
V2 = v(R) for the red, together with B\ = Bu, B2 = Br. 
With the mean wavelengths of t he filters, 2942.8 A for 
F300W and 6162. 8A for F622W (|Biretta et al. 2000D , a 
cut on the signal-to- noise ratio S/N > S/N cut of the input 
images limits the spectral index error to at most 



< 



1 



In 



6162.8 
2942.8 



1.9 



S/N min S/N n 



(3) 
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Fig. 1. The jet in red light (620 nm) 
after background subtraction. 
Logarithmic grey-levels run from 
to 0.04/xJy/pixel, 0'.'08 effective beam 
size, CK'045 pixel size. The quasar 
core lies 10" to the northeast from 
A. The labelli ng of the jet featur es as 
introduced by Lclicvrc ct al. (1984) an d 
extend ed by Roser & Meisenheimei 
(1991), together with the hot spot 



nomen clature from Flatters & Conway 
(1985| ) is also shown. Note that the 
labelling used by Bahcall et al. (1995) 
is slightly different. 



Fig. 2. The jet in UV light (300 nm) 
after background subtraction. 
Logarithmic grey-levels run from 
to 0.014/xJy/pixel, 0"06 effective 
beam size, 0'.'045 pixel size. 



In our case, the UV-band S/N is always much inferior to 
that in the red band, so the error will actually be domi- 
nated by the UV noise and hence smaller than the max- 
imum from Eqn. ^ for the most part. The noise in the 
input images is calculated from the calibrated images and 
includes shot noise due to both observed photons and dark 
current as well as the read noise, but excludes systematic 
errors, whose expected magnitude we now examine. 



5.2. Systematic errors for a spectral index 
determination 

We consider the systematic errors which may be intro- 
duced when combining two images taken through differ- 
ent filters, or by different instruments and telescopes. The 
main danger in a determination of spectral gradients lies in 



a misalignment between the images which would introduce 
spurious gradients; referring the flux to different effective 
beam sizes will lead to wrong spectral index determina- 



tions as well. Following considerations given in Sect. A.l 



we deduce that a 5% limit on the flux error due to mis- 
alignment requires aligning the images to better than 10% 
of the effective PSF full width. The error in the PSF deter- 
mination is negligible when, as in our case, the smoothing 
Gaussian is much wider than the PSF. 

The 0'.'2 effective resolution aimed for thus requires 
knowing the relative alignment of all images in the data 
set to better than 20 mas or 0.44 PC pixels. The abso- 
lute telescope pointing does not need to be known for this 
purpose as we tie all positions to the quasar core as origin. 

After detailed investigations of all issues rela ted to rel- 
ative alignment of HST images (see Sect. A. 2 ), we used 
the following procedure: as the first step, all exposures 
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-1.5 
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Fig. 3. Optical spectral index anu 
(f(u) oc v a ) at 0'.'2 resolution. 
The contours show the red-band im- 
age and are logarithmic with a fac- 
tor V2 from 30 /^Jy/beam to 340 
/iJy/beam. 



through one filter and within one set of observations are 
summed up using the relative positions from the "jitter 
files" provided as part of the observing data package (typ- 
ical error 5 mas). The exposures through F300W were 
distributed over three different sets of observations, so we 
obtain three intermediate images. For aligning these with 
each other, we use the average of the positional shifts of 
the quasar and three astrometric stars (observed on the 
Wide Field chips) measured on the intermediate image 
or the short exposure in each set of observations (error 
10 mas). The field flattener windows inside WFPC2 intro- 
duce a wavelength dependence of the plate scale, which 
has to be removed prior to combination of the images 
taken through different filters. We therefore sample all 
four intermediate images (three UV and one red) onto a 
grid of one tenth of their average plate scale with a pixel 
size of 0'.'0045548. We sum the UV intermediate images 
to give a UV final image. Because of the scale change, we 
have to align this to the red image using the quasar posi- 
tion only (estimated error 10 mas-15 mas). We then rebin 
both images to the final pixel size of 0'.'045548. Adding up 
all errors in quadrature, the error margin of 20 mas is just 
kept. 



Note that the rebinning has only matched the scales of 
the two images, but not removed the geometric distortion 
of the focal plane. The next-order wavelength-dependent 
term in the geometric distortion solution is at least two or- 
ders of magnitudes smaller than the scale difference. Both 
the red and the ultraviolet image can thus be assumed to 
have identical distortion solutions. We therefore ignore the 
geometric distortion for the remainder of this paper, but 
caution that it must be taken into account when compar- 
ing these data with data obtained at other instruments. 



5.3. Calculation of the optical spectral index map 

We have now aligned the images in the two filters to better 
than 0'.'02. When comparing flux measurements on these 
images, they have to be referred to a common beam size. 
The maps are therefore smoothed to a common 0'.'2 using 
Gaussians with a width matched to the effective FWHM of 
the PSF in each image. This operation retains the original 
pixel size. As few point sources are available, the effective 
FWHM is determined on PSF models generated using the 



TinyTIM software flKrist 1999| ), resulting in 0'.'08 on the 
red and 0'.'06 on the UV-band image. Because of sampling 
effects, the measured widths are larger than the size of 
the Airy disk of the HST's primary at the corresponding 
wavelengths. 

Fig. H shows the spectral index map obtained from the 
images (Figs. [j] and ||) according to Eqn. [j]. Only those 
points are shown in the map which have an aperture S/N 
of at least 5 on both images. This limits the error in the 
spectral index to 0.4 (Eqn. |J). The error is less than 0.1 in- 
side all knot regions, so that a colour change indicates sta- 
tistically significant variations of the spectral index. Fig. |^ 
shows the spectral index along a tracing of the jet obtained 
from a rotated map with the jet's mean position angle of 
222?2 along the horizontal. 



5.4. Colours of the extensions 

The outer extension is not bright enough to show up on 
the map in its entirety. The variation of olru across it 
are consistent with a spiral galaxy. The inner extension's 
two knots have markedly different colours on the spectral 
index map. Knot In2 shows a spectral index gradient, a sa 
—0.9 to —0.3 roughly parallel to the jet and outwards from 
the quasar position. Knot Inl has a much steeper spectral 
index of about —1.5 and shows no gradient. Interestingly, 
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Fig. 4. Run of the red-band bright- 
ness and optical spectral index along 
the outer half of the jet in 3C273, 
for a 0'.'2 beam, olru was deter- 
mined from Figures 1 and 2, while 
ctBRi for a l'/3 beam is taken from 
Roser & Meisenheimer (1991). For 



comparison, we show a&t, the cor- 
responding spectral index obtained 
from synchrotron spectra fitt ed by 
Meisenheimer ct al. (1996a| ) and 
Roser et al. (2000|) . While the ob- 



served spectral index agrees with 
older data, it is now clear that the fit 
is inadequate for the optical part of 
the spectrum. The steeper spectral 
index of the fit may either be due 
to contamination of the infrared flux 
by a "backflow" component around 
the jet, or because of the presence of 
a second high-energy electron pop- 
ulation in the jet which is not in- 
cluded in the fit (see Sect. 6.3. l| for 
a discussion of the discrepancies) . 



the southern extension S has the flattest spectrum (olru 
—0.4) of all regions on the map. 



5.5. Spectral index of the jet 

The optical spectral index declines globally outwards from 
—0.5 near the onset of the optical jet at A to —1.6 in 
D2/H3. This trend does not continue into the hot spot, 
again stressing the physical distinction between jet and 
hot spot. The general steepening is in agreement with pre- 
vious determinations of the knots' synchrotron spectrum 
which showed a decline of the cutoff frequen cy outwards 
( Meisenheimer et al. 1996a ; Roser et al. 2000 ) . The optical 
spectral index olbri determined at l'/3 resolution (Roser 
& Meisenheimer 1991) agrees very well with our new de- 
termination of olru at the much higher resolution (Fig. ||). 
The only dis crepa ncy arises in C2/D1; we defer a discus- 
sion to Sect. 3.3.1. The smooth variations of the spectral 



index along the jet show that the physical conditions in 
the jet change remarkably smoothly over scales of many 
kpc. There is no strict correlation between red-band sur- 
face brightness and spectral index like that found in the jet 
in M87 ( [Meisenheimer et al. 19964 |Perlman et al. 200 ip . 

There is a marginally significant flattening of the spec- 
trum in the transitions A-Bl, B1-B2, C1-C2, and moving 
out of C2, consistent with the absence of losses. In any 
case, the overall steepening of the spectrum (from region A 
down do D2/H3) is less rapid than that within individual 
regions (e. g., A and B2). There is no significant steepen- 
ing from Dl out to the bridge between H3 and H2, despite 
large variations in surface brightness. 



The criss-cross morphology in CI, C2, Bl, B2 and Dl 
is reflected on the spectral index map as a band of one 
colour crossing a second one. This is most clearly seen in 
knot CI which has a green band of a sa —1.1 across an 
orange region of a sa —1.4, supporting the interpretation 
of two emission regions appearing on top of each other. 
The spectral index near the hot spot shows a flip from 
flat (-1.0) to steep (-1.5). 

6. Discussion 

6.1. The jet's extensions 

Before analysing the jet emission proper, we consider the 
relation of the extensions to the jet. Our new images con- 
firm the outer extension's nature as spiral galaxy. The 
spectral index map gives us few hints about the nature of 
the emission mechanism for the inner and southern ex- 
tensions. We have therefore measured their integral fluxes 
at 300 nm and 620 nm (this paper) and at 1.6/x (HST 
NICMOS camera 2 imaging^, in prep.). Fig. [| shows the 
resulting spectral energy distributions. Inl, with the steep- 
est olru, is also the brightest in the infrared. The infrared 
flux points of both extension S and of In2 lie on the con- 
tinuation of the power law with index —0.4, determined 
on the optical spectral index map. If this power law were 
valid up to radio wavelengths, it would predict a flux at 
A6cm of about 50/iJy for the northwestern knot Inl, and 
of about 45/i Jy for the southern extension S. This is a fac- 
tor of order 100 below the flux from the jet at this wave- 



length (Conway et al. 1993), and any emission at this level 
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would not be detected even on our new VLA maps with 
an RMS noise of » 0.8 mJy. There is no similarly obvious 
extrapolation of Inl's SED to radio wavelengths. A full 
comparison will be presented in a future paper. 

The inner extension shows a polarisation signal on un- 
published ground-based maps made by our group and 
on polarimetric images obtained with the Faint Object 
Camera on board HST (Thomson ct al. 1993). A polarisa- 



tion signal could be due to scattered quasar light ( Rbser & 



Meisen heimer 1991) or, together with possible radio emis- 
sion, synchrotron radiation. We aim to clarify the nature 
of the inner extensions by planned polarimetric and spec- 
troscopic observations with the VLT. 

6.2. The hot spot 

The optical counterpart to the radio hot spot (H2) ap- 
pears very faint on both images, and the spectral index 
map shows a flip fro m flat to steep there. Both feat ures 
are explained by the Meisenhcimcr fc Heavens (1986 ) hot 
spot model: the lifetime of electrons emitting in the op- 
tical is quenched by the strong magnetic field in the hot 
spot. The spectral index flip is expected if there is an off- 
set between the emission peaks at different frequencies. 



This offset is predicted by the theory by Meisenheimer & 
Heavens (1986| ) and is clearly seen when comparing opti- 



cal and radio images (Roser et al. (1997); a detailed com- 
parison of all data will be presented in a future paper). 
By its spectrum and morphology, the hot spot appears 
exactly as expected for a high-loss synchrotron emission 
region downstream of a localized strong shock with first- 
order Fermi acceleration. Finally, we note that although 
the term "hot spot" is historically in use for this part of 
the jet, the term is ill-defined and may have become in- 
adequate, but we defer a discussion of its adequacy to a 
future paper. 

6.3. Main body of the jet 

Our HST images (Figs. [J and |J) show a close coincidence 
of the jet's morphology over a factor of 2 in frequency 
(optical and UV). The spectral index map (Fig. ||) shows 
amazingly smooth variations in the physical conditions 
over the entire jet. The electrons' synchrotron cooling 
leaves little imprint on the spectral index along the jet, 
contrary to expectations. 

6.3.1. The spectral index at HST resolution 

The two striking features of our spectral index map are 
the smooth variation of the spectral index over the entire 
jet (Fig. |J), and the lack of a strong correlation between 
red-band brightness and sp ectral index (Fig. |g|) like th at 
observed in the jet o f M87 ( Meisenhcimcr ct al. 1996a , g; 
Pcrlman ct al. 2001 ). There are large but smooth varia- 



tions of the spectral index, while the surface brightness 
remains fairly constant over the jet's projected extent of 
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Fig. 5. Broad-band flux points for the radio-quiet exten- 
sions to the jet from HST imaging at 1.6 /i, 620 nm and 
300 nm. 
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Fig. 6. Optical spectral index against red-band brightness 
for the outer half of the jet in 3C 273 (data as in Fig. [|). 
There is no strict correlation like that observed in the jet 
of M87. 

about 10". Conversely, there are large local variations of 
surface brightness without strong changes in the spectral 
index. We do not observe abrupt changes of the spectral 
index inside the jet like we do in the hot spot. The run 
of the spectral index is thus consistent with the complete 
absence of energy losses over scales of many kiloparsecs, in 
spite of the large observed synchrotron luminosity, which 
indicates the need for re-acceleration of particles inside 
the jet. The mismatch between synchrotron cooling scale 



and extent of the jet even led Guthrie & Napier (1975) to 
conclude that the jet's emission could not be synchrotron 
radiation. 

The same lack of apparent synchrotron cooling is im- 
plied by the surprising overall correspondence between 
o-bri and am/ (Fig. |4|) which have been determined at 
vastly different resolutions (l'/3 and 0'.'2, corresponding 
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to 4.2/i6okpc and 640ft.6oPCj respectively). The small dif- 
ferences (strongest in regions A and C2/D1) can be ex- 
plained at least in part by the different beam sizes: the 
spectral index determined by smoothing the HST im- 
ages to l'/3 makes the discrepancy smaller. In addition, 
the emission region becomes wider towards longer wave- 
lengths, an effect already noticeable in the near-infrared 
K' band (Neumann et al. 1997) and attributed to a "back- 
flow" (Roser et al. 1996). The presence of such a steep- 
spectrum, diffuse component more extended than the jet 
channel might also explain the discrepancy. 

The general agreement indicates a constancy of the 
shape of the spectrum even at the scales resolved by HST. 
To first order, the appearance of the jet at all wave- 
lengths can be explained through such a constant spec- 
tral shape consisting of a low-frequency power law and a 
high-frequency curved cutoff. The assumed spectral shape 
is shifted systematically in the log v — log S v plane, pass- 
ing different parts of the curved cutoff through the optical 
flux point. This leads to similar morphological features 



at all wavelengths, in accord with observations (Conway 



et al. 1 393; Bahcall ct al. 1995). Moving the spectral shape 
through a constant optical flux point leads to a correla- 
tion of steeper optical spectral index with higher radio 
flux. This correlation reproduces the observed tenfold in- 
crease of the radio surface brightness from A to D2/H3 
( |Conway et al. 1993 ) and the overall steepening of the op- 
tical spectral index. However, this steepening happens on 
much longer time scales than expected from synchrotron 
cooling alone. 

Assuming an electron energy distribution with a fairly 
sharp cutoff, drastic jumps in ajm would be expected at 
the locations of the shock fronts if the jet knots were 
(strong) shocks like the hot spot without extended re- 
acceleration acting between them. The absence of strong 
cooling makes it impossible to pinpoint localised sites at 
which particles are either exclusively accelerated or exclu- 
sively undergo strong losses. Any acceleration site must 
therefore be considerably smaller than the beam size we 
used (640/i^q 1 pc) , and these acceleration sites must be dis- 
tributed over the entire jet to explain the absence of cool- 
ing. This, together with the low-frequency spectral index 
of « —0.8, corresp onds to the jet-like acceler ation mecha- 
nism proposed by Meiscnhcimcr ct al. (1997). 

Meisenheimer ct al. (1996a) have presented fits of syn- 
chrotron continua to the observed SEDs of knots A, B, C, 
D and the hot spot H at l'/3 which have also been used 
by Roser et al. (2000). It is noted that the optical spec- 
tral index predicted from these spectra (afn in Fig. |^) is 
always steeper than the observed spectral index, that is 
below asm and anu- This indicates that the fitted spec- 
trum is not fully adequate at the highest frequencies. This 
may again be due to contamination of the near-infrared 
flux by the same "backflow" component mentioned above. 
The contamination would make the IR-optical spectral in- 
dex (which dominates the run of the spectrum at high 
frequencies) steeper than the optical spectral index, as 
observed. Alternatively or additionally, there could be de- 



viations of the optical-UV spectral shape from a standard 
synchrotron cutoff spectrum which may be interpreted 
as the first observational hint towards the existence of a 
second, higher-energy electron population producing the 
X-ray emission (Roser et al. 2000). Detailed statements 
about the adequacy of model spectra have to be deferred 
to a future paper considering the full radio, infrared and 
optical data set. 

6.3.2. Can beaming account for the lack of cooling? 

Heinz & Bcgelman (1997") proposed that sub-equipartition 
magnetic fields combined with mildly relativistic beam- 
ing could explain the lack of cooling in the jet of M87 - 
which is, however, ten times shorter than that of 3C273. 
As an alternative to postulating re-acceleration, we con- 
sider whether low magnetic field values and beaming could 
lead to electron lifetimes sufficient to allow electrons to be 
accelerated at region A to illuminate the entire jet down 
to the hot spot over a projected extent of 32/i f r 1 kpc (the 
argument will become even more stringent by demanding 
acceleration in the quasar core). We consider the electron 
lifetime against synchrotron and inverse Compton cool- 
ing off cosmic microwave background photons; the syn- 
chrotron s clf-Compton proc ess is negligible for electrons 
in the jet ( Roser et al. 2000 ), as is Compton scattering off 
the host galaxy's star light. 

The total energy loss rate of an electron with energy E 
due to synchrotron radiation and inverse Compton scat- 
tering, averaged over many pitch-angle scattering events 
during its lifetime, is 

2 

(4) 



dE 
dt 
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where [/tot = Ucbk{ z ) + C/ m ag is the sum of the energy 
densities of the background radiation and magne tic field, 
respectively, and ctt is the Thomson cross-section ( Longaii 
1994 ). We integrate this equation from E = oo at t = 
to E(t), assuming (3 = 1 (appropriate for the highly rel- 
ativistic electrons required for optical synchrotron radia- 
tion) and substitute for the electron's energy E = jm e c 2 . 
Inverting yields the maximum time that can have elapsed 
since an electron was accelerated, given its Lorentz factor 
7 (van der Laan & Perola 1969): 



(5) 



3 -crTc£/tot7 

This is the "electron lifetime", inversely proportional to 
both the energy density in which the electron has been 
"ageing" , and the electron's own energy. 

As most of the electron's energy is radiated at the syn- 
chrotron characteristic frequency oc 7 2 -B, we can substi- 
tute for 7 in Eqn. in terms of the observing frequency 
and the magnetic field in the source. Hence, Eqn. ^| be- 
comes (in convenient units) 
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where -B_g,j C t is the magnetic flux density in nT of the 
jet field, the background radiation energy density has 
been expressed in terms of an equivalent magnetic field, 
B-g ic = (1 + z) 2 x 0.45 nT, and ^15 is the observing fre- 
quency in 10 15 Hz (van der Laan & Perola 1969). Note 
that as the substituted 7 oc B~?, setting B_g,j C t = is 
now meaningless. 

To be fully adequate for electrons in a relativistic jet 
at cosmological distances, the equation needs to be mod- 
ified. Firstly, the frequency local to the source is (1 + z) 
times the observing frequency because of the cosmologi- 
cal redshift. Furthermore, the radiating electron may be 
embedded in a relativistic flow with bulk Lorentz fac- 
tor r with three consequences: relativistic time dilation 
and Doppler shift, and boost of the background radia- 
tion energy density. The relativistic time dilation enhances 
the electron lifetime in the jet frame by a factor T. The 
Doppler shift between the emission frequency Vi ni in the 
jet frame (equal to the characteristic frequency) and the 
observation frequency i/ b s is given by is- m t = v^U^ 1 , 
where V(T,9) = [T(l - (3- iCt cosfl)] -1 , the Dop pler boost- 
ing fac tor fo r an angle 9 to the line of sight (see Hughes fc 
e. g.). A relativistic flow perceives the energy 
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density of the backgroun d radiation field boosted up by a 
factor T 2 (Dcrmer 1995). Inserting these yields 
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The electron lifetime attains a maximum value at a 
certai n value of the j et's magnetic field (Fig. |t|) ( |van der 
Laan fc Perola 1969). On either side of the maximum, 
the lifetime is decreased by larger losses suffered before 
it is observed: for a higher magnetic field in the jet, the 
synchrotron cooling is more rapid (t cx 1/iTtot, Eqn. |5|). A 
lower jet field requires an electron of higher Lorentz factor 
for emission at the given frequency, which also suffers more 
rapid losses (t oc I/7, Eqn. ||) . By differentiation of Eqn. fj], 
the value of the maximum lifetime is 



t-maxi.Bic'D, z) 
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Note that the largest possible value for the factor \jT>jY 
is y/2. imax is a firm upper limit for the lifetime of a 
synchrotron-radiating electron from a source at redshift 
z in a flow of bulk Lorentz factor T, whatever the mag- 
netic field strength in the source. It is deduced only from 
the fact that synchrotron emission is observed at a certain 
frequency, and that electrons which can radiate at this fre- 
quency suffer drastic energy losses either by synchrotron 
or by inverse Compton cooling between the time of ac- 
celeration and the time of emission. The only additional 
assumptions are rapid pitch-angle scattering and a homo- 
geneity of conditions throughout the electron's lifetime. 

The VLBI jet close to the core has a line-o f-sight angle 
near 1 0° and a bulk Lorentz factor near 10 ( Abraham fc 
Romem 1999). A line-of-sight angle 9 ss 45° has been 
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Fig. 7. Maximum age against synchrotron and inverse 
Compton cooling (off microwave background photons at 
3C 273's redshift) of an electron radiating at observed UV 
wavelengths plotted against the jet magnetic field. Solid 
line, no beaming; dashed line, relativistic Doppler boost- 
ing with r = 10 and a line-of-sight angle 9 = 45°. 



inferred for the flow into the hot spot from indepen- 
dent considerations of the jet's p olarisation change there 
and the hot spot's morp hology ( Conway fc Davis 1994 ; 
[Mciscnhcimer et al. 1997| ) . We have plotted the lifetime of 
an electron responsible for emission from the jet in 3C273 
observed at 300 nm as function of the jet's magnetic field 
in Fig. ^ for the extreme cases of no beaming in the op- 
tical jet and beaming identical to that in the VLBI jet 
with r = 10, 6 = 10°, and for an intermediate case with 
T = 10,6 = 45° (though note that T — 10,9 — 10° is un- 
realistic as there is a difference in position angle between 
the VLBI jet at 244° and the arcsecond jet at 222°). 

The equipartition flux densities derived for the jet lie 
in the range of 15 nT ( Neumann 1995| ) up to 67 nT for 
region A (Roser ct al. 2000), leading to maximum ages of 
100-800 years, less than the light travel time from one 
bright region to the next. In the absence of beaming ef- 
fects, the largest possible lifetime for electrons in 3C 273 
from Eqn. || is about 58 000 y, again short of the required 
values. The "boosted lifetime" can be at most V2 larger 
than this. There is thus no combination of T, 9 which en- 
hances the electron lifetime to the 100 000 y required for 
illumination of the entire jet in 3C 273 by UV-radiating 
electrons. 

Thus, the invocation of mild or even drastic beaming 
and/or sub-equipartition fields cannot resolve the discrep- 
ancy between the synchrotron loss scale and the extent of 
the o ptical jet of 3C 273, as has been possible for the jet in 
M87 ( Heinz fc Begclman 1997 ). As another alternative to 
invoking quasi-permanent re-acceleration, the existence of 
a "loss-free channel" in which electrons can travel down 
a jet without synchrotron cooling has been proposed by 
|Owcn ct al. (1989 ). As an extreme version of this case, 
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we assume that the electron travels along the jet in zero 
magnetic field and is observed as soon as it enters a fil- 
ament with magnetic field Bfn. The energy loss between 
acceleration and synchrotron emission is then only due to 



the hot spot. Recently, the detection of extended X-ray 
emission from the jet in PKS 0637—752 by the new X- 
ray observatory Chandra has been reported and a simi- 



lar explanation has been brought forward (Chartas et al 



Compton scattering. The lifetime in 3C 273 is then 2000; 


Schwartz et al. 2000; 


Tavccchio et al. 200C : 


Celotti 


y xr~ i v /B_9, fi iP/i/i 5 . Again, if the jet flow in et al. 2001 


). Like for 3C273, in situ re-acceleration is re- 



3C 273 is highly relativistic, the electrons suffer heavy in- 
verse Compton losses and the lifetime mismatch persists. 
In any case, it remains to be shown that the "loss-free 
channel" is a physically feasible configuration of an MHD 
jet. 

7. Conclusion 

We have analysed deep HST images at 300 nm and 600 nm 
(Figs. and ||) of the jet in 3C 273 and constructed an op- 
tical spectral index map at / .'2 resolution (Fig. |||). The 
optical spectral index varies smoothly over the entire jet, 
indicating a smooth variation of the physical conditions 
across the jet. Unlike in M87 ( Meiscnhcimcr et al. 1996a; 
Pcrlman et al. 2001), there is no strong correlation be- 
tween optical brightness and spectral index (Figs. [| and 
^) . The spectral index map thus shows no signs of strong 
synchrotron cooling at any location in the jet. Particle ac- 
celeration at a few localised sites in the jet is not sufficient 
to explain the absence of strong cooling. This does not pre- 
clude the possibility that the enhanced-brightness regions 
are shocks - but even if they are, re-acceleration between 
them is necessary to explain the observed spectral index 
features. We have further shown that relativistic effects 
cannot lead to significant enhancements of the electron 
lifetime in 3C 2 73's jet, whatever the bulk Lorentz fac- 
tor (Sect. |6.3.2 ), strengthening previous electron lifetime 
arguments. The need for a continuous re-acceleration of 
electrons emitting high-frequency synchrotron radiation 
in the jet of 3C 273 is thus evident. 

Mechanisms have been proposed which can explain the 
apparent lack of cooling by distributed re-acceleration. 
These i nclude acceleration by reconnection in thin fila- 
ments ( Lesch fc Birk 1998 ) and turbulent acceleration 



( Manolakou et al. 1999 ). Both processes manage to main- 
tain the injection spectrum over distances much larger 
than the loss scales, although the latter so far only main- 
tains cutoff frequencies in the range of 10 12 Hz-10 13 Hz, 
i. e., below the values observed in 3C273. 

We note that for those jets which are bulk relativis- 
tic flows at high Lorentz factors, the increased inverse 
Compton losses form a further sink of energy that has 
to be filled by re-energization processes inside the jets. 
This requirement becomes more severe at higher redshifts. 



As has been suggested previously (Celotti et al. 2001) 
inverse Compton scattering off cosmic microwave back 
ground photons might explai n the so far unaccounted-for 



X-ray flux from 3C 273's jet flRoscr et al. 2000 ; [Marshall 



ct al. 2001 ) and from Pictor A's jet and hot spot ( [Wilson 
ct al. 2001 ). If this is true, the outward-decreasing X-ray 
flux from the jet in 3C 273 indicates that the jet is still 
highly relativistic near region A and slows down towards 



quired in PKS 0637—752 to explain the mismatch between 
de-projected extent of the jet (> 1 Mpc) and inverse- 
Compton loss scale (lOkpc) (Tavecchio et al. 2000). 

In order to understand the physical conditions in ex- 
tragalactic jets like that in 3C 273, it is necessary to detect 
even the most subtle variations in the parameters describ- 
ing the synchrotron spectrum (that is, cutoff frequency, 
break frequency and especially low-frequency spectral in- 
dex), which requires the deepest images at the highest 
resolutions and in many wavelength bands to detect vari- 
ations at all. A theory of the physical processes at work 
in this jet initially has to explain both the simple overall 
spectral shape, as well as its constancy over scales of many 
kiloparscc. The details of this physical process will be con- 
strained by the subtle deviations from the simple spectral 



shape, such as those tentatively identified in Sect. 6.3.1 
We aim to find these deviations with the full data set, in- 
cluding new radio and near-infrared data in addition to 
the presented optical and UV images. 

Acknowledgements. We are grateful to D. van Orsow for his 
assistance with the HST observations. This research has made 
use of NASA's Astrophysics Data System Abstract Service. 

Appendix A: Image alignment 

A.l. Flux errors from misalignment 

Consider two images which are not registered correctly 
and which have slightly different PSF widths. When per- 
forming photometry on these images, we assume that they 
are registered perfectly and have identical, known beam 
sizes. This amounts to making a flux measurement in a cer- 
tain aperture in one image, but in a slightly offset aperture 
of slightly different size in the second image. The error is 
largest in the steepest gradients in the image, which are 
the flanks of the PSF of width a, smoothed to the desired 
effective beam size a e ff, with a 2 s = c 2 +Cg mooth . Its magni- 
tude can be assessed by considering the PSF as a Gaussian 
at given position and of given width, and as smoothing fil- 
ter a second Gaussian slightly displaced from the PSF and 
of width slightly different from that achieving the desired 
effective beam area. The result of the performed wrong 
flux measurement is then proportional to the integral 



exp 



x 2 + y 2 
lia 2 

zo PSF 



exp 



(x — Sx) 2 + y 2 
2(cr smoo th + So-) 2 



dxdy, 



where Sx is the offset of the aperture from the correct 
position and 8a is the error in the determination of the 
PSF width. The correct measurement is obtained by set- 
ting 8x — and 5a — 0, and from this one obtains the 
fractional flux error as a function of the two errors. The 



S. Jester et al.: HST optical spectral index map of the jet of 3C273 



11 




Fig. A.l. The separation of a number of stars was mea- 
sured both on a F336W (UV) and a F673W (red) im- 
age of R136. This plot shows the difference in separation 
between the frames, AS, plotted against the separation 
on the "red" image, S. The best-fitting straight line is 
shown. AS grows systematically with S, indicating a dif- 
fering scale between the frames. The slope of the line is 
(1.14 ± .04) x 10~ 3 , the intercept with the ordinate is 
(0.3±8) x 1CP 2 . The slope is not changed within the errors 
by forcing the line to pass through the origin. 



fractional flux error A/ for a misalignment is 



Af = l-e 2 ° 
Sx 2 



^ ^smooth 
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Hence, the relative flux error is better than 5% if the mis- 
alignment Sx < 10%. Similarly, for a wrong PSF width, 
the fractional error is 



2S(J ° 2 •* x 

A J = — ll OCT < CT smoot h- 

^smooth ^ c ff 



(A.3) 



This error is negligible if the desired effective PSF is much 
larger than the intrinsic PSF of the input images, as is the 
case here. 



A.2. Refined image alignment 

The two obvious ways to determine the relative shift be- 
tween any two images are measuring the positions of point 
sources on the various images, or using the engineering 
files (also termed jitter files) provided as part of the ob- 
serving data package. If the telescope pointing was known 
to be precise to better than 0'.'02, we could simply rely on 
the commanded shifts. 



The HST has a pointing repeatability within a single 
telescope visi^ of about 5 milli-arcsecond (mas). The off- 
sets between individual exposures are accurate to 15 mas, 
leading to a total error of about 16 mas. This is already 
comparable to the demanded accuracy. The situation is 
expected to be worse for relating exposures in different 
visits, when the telescope has been pointing elsewhere in 
the mean time. 

Unfortunately, measuring the positions of only a few 
(four, in our case) astrometric reference stars does not 
immediately lead to accurate measurements of the tele- 
scope's pointing - especially since in our case, each star 
lies on a different chip. The undersampling of the telescope 
PSF by the WFPC2 pixels and the so-called sub-pixel 
scattering of the WFPC2 detectors lead to an additional 
scatter of the centroid positions of point sources, approx- 
imately uniformly distributed between +0.25 and —0.25 
pixels and in excess of statistical uncertainties (Lallo, 
1998, priv. comm.). With only a small number of centroid- 
able point sources available, the centroiding errors are of 
the same magnitude as the intrinsic pointing errors of the 
telescope. 

It is therefore worth considering the pointing error 
sources en detail to ensure that the alignment is at the 
required 0'.'02 level^j. Alignment errors can be caused by 
roll or pointing errors and less obviously by a scale differ- 
ence between exposures using different filters in the same 
camera. The importance of the various alignment error 
sources can be estimated by considering the effect they 
have on the hot spot location if the quasar images are as- 
sumed to coincide. The hot spot is separated by 22" from 
the quasar, corresponding to about 480 PC pixels. 

The roll repeatability of the HST is at the 10" level. 
The engineering data provided with HST exposures record 
roll angle differences of about 6". This is well below the 
rotation of 3.'5 which would produce a 0.44 pixel difference 
over 480 pixels. The telescope roll differences can thus be 
neglected. 

The engineering files record the telescope pointing in 
three-second intervals and can be used to calculate the 
offsets. Their accuracy is only limited by the so-called 
"jitter" , vibrations due to thermal effects. The jitter was 
below 10 mas in all exposures, and below 5 mas in most. 
There is an additional uncertainty from the transforma- 
tion between the telescope's focal plane and the detector: 
the location of a camera inside the telescope may change 
slightly over time (shifts, rotations, or both). This uncer- 
tainty is irrelevant for relative positions as long as the lo- 
cation and orientation of WFPC2 and the Fine Guidance 
Sensors (FGS, these perform the guiding observations) in- 



5 A visit is a set of exposures that should be observed to- 
gether; the HST equivalent of an observing run. The telescope 
is aligned to its guide stars only at the beginnin g of a visit. 

6 Th e use of the quasar's diffraction spike by Bahcall et al 



(1995 ) to locate the QSO's center is only expected to work 
well for sources near the camera center (Krist, priv. comm.). 
We failed to reproduce their quoted accuracy with our data. 
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side the telescope is stable, which is the case for the em- 
ployed shifts of about 1" and for the timescales between 
the visits. 

Within a single visit, the engineering file information 
is used to obtain relative offsets, with a typical 5 mas 
error. The values differ from the commanded shifts by a 
few milli-arcseconds at most. All of these offsets are by an 
integer number of PC pixels. 

Between the various visits, the telescope has been 
pointing to a different part of the sky. One should there- 
fore not assume that the relative shifts between various 
visits as determined from the jitter files are as accurate 
as shifts within one visit. We therefore measured the po- 
sitions of the four astrometric reference stars on each of 
the short exposures (in fact, one of the stars is very faint 
in the UV, so the position was determined on a long ex- 
posure for this one). The shifts determined from the four 
point sources' positions have a typical standard deviation 
(accuracy) of 15 mas in each coordinate. On the three vis- 
its' sum images aligned this way, the scatter of the quasar 
image position is less than 5 mas and 7 mas in x and y, re- 
spectively. This means that although the measured shifts 
have a fairly large scatter, the resulting value is precise to 
about 10 mas. Finally, we note that the observed shifts 
of the stellar positions and those obtained from the jitter 
files agree to better than 0'.'02 in all cases, with an RMS 
value of O'.'Ol. There are, however, systematic differences 
between these and the commanded values. Hence, we do 
not blindly rely on the latter. 

Because of differential refraction i n the MgF2 field flat - 
tener windows employed in WFPC2 ( Frauger et al. 1995| ), 
the pixel scales of images taken through the F622W and 
F300W filters differ by about 0.1%. This alone is enough 
to eat up the alignment error budget of 0.44 pixels over 
480 pixels separation. The wavelength dependence was ex- 
pected from ray-tracing studies of the WFPC2 optics. Its 
presence and magnitude were experimentally confirmed by 
comparing archival images^] of the star cluster R136 taken 
through similar filters as those employed in the present 
work (F336W and F672N) (Fig.[0]). The scale difference 
has to be removed before combining the images to a spec- 
tral index map. This was done in the following manner: the 
plate scale of each image was calculated using the param- 
eters in Frauger et al. (1995| ). All images were resampled 
to a grid with pixel size of 0'.'0045548, which is one tenth 
of the average of the original scales, using bilinear inter- 
polation. The result was then binned in blocks of 10 x 10 
pixels to a common pixel size of 0'.'045548. The scale of 
the two images is then identical to better than 1 part in 
10 000. 

The resampling required the use of the quasar image 
as common reference point between the two filters. The 
QSO's position can be determined to about 10-15 mas by 
centroiding routines on the unsaturated images. 



7 Data from proposal 5589, PI John Trauger 
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